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Abstract 



We study the utility of a complex Langevin (CL) equation as an alternative for the Monte 
Carlo (MC) procedure in the evaluation of expectation values occurring in fermionic many- 
body problems. We find that a CL approach is natural in cases where non-positive definite 
►^ ■ probability measures occur, and remains accurate even when the corresponding MC calculation 

develops a severe "sign problem" . While the convergence of CL averages cannot be guaranteed 
fSJ . in principle, we show how convergent results can be obtained in three examples ranging from 

f^ ' simple one-dimensional integrals over quantum mechanical models to a schematic shell model 

0^ ' path integral. 
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1 Introduction 

There has been recent significant progress in the large-scale numerical computation of nuclear prop- 
erties in the shell model, using the Auxiliary Field Path- Integral (AFPI) Monte Carlo method ||I|. 
Due to the benign scaling of the computational effort with the single-particle basis, properties of 
large nuclei can now be calculated that are out of the reach of conventional diagonalization methods 

i|. 

Quite generally, AFPI treatments of interacting fermion systems with the Monte Carlo (MC) 
method are difficult for certain realistic Hamiltonians, and for odd-particle configuration^ This 
difficulty, also known as the "sign problem" , is the prime impediment to large scale computational 
efforts both in nuclear and condensed matter physics calculations. Briefly, repulsive interactions 
and/or odd-particle configurations can lead to probability distributions (integration measures for 
the auxiliary fields) that are negative, or even complex. As the MC update algorithm relies on a 
positive-definite measure, the sign of the distribution is made part of the observable being calculated. 
Under circumstances where the average sign of the distribution is small, the expectation value is the 
ratio of two very small numbers, that converges only asymptotically. 

A number of alternatives to or extensions of the MC method have been proposed over the years, 
among them hybrid methods combining the MC method with replication mechanisms for importance 
sampling H, random- walk branching H , and diagonalization over optimal bases using variational 
techniques^]. Here, we investigate the possibility of replacing the MC method altogether with one 
based on the complex Langevin (CL) equation, at least in those cases where the sign problem is 
prominent. The CL equation has received considerable attention in connection with lattice gauge 
theory calculations, where either static charges or a non-zero chemical potential give rise to complex 
actions pi . It has been abandoned mostly because of the perception that complex Langevin averages 
ought not to be trusted, due to the possibility that they can be non-stationary. Here, we show that 
the problem of non-stationarity can often be avoided, and does not play a role in a large class of 
Hamiltonians which give rise to Langevin equations with fixed points in the complex plane. 

By way of introduction and to establish notation, we briefly review the AFPI Monte Carlo method 
in the next section. Section 3 discusses the application of the CL method to simple one-dimensional 
integrals abstracted from those occurring in fermionic many-body systems, while Section 4 expands 
this to a simple toy Hamiltonian with characteristics reminiscent of shell models. In Section 5 we 
apply the method to the Lipkin model (the MC treatment of which is similar in character to the full 
shell model), and close with conclusions, in Section 6. 

2 Auxiliary Field Path-Integral Monte Carlo 

In the AFPI method, the significant savings in computational effort are obtained through a Hubbard- 
Stratonovich (HS) transformation [B| (see below) of the imaginary-time evolution operator 

t/ = exp(-/3iJ) . (2.1) 

The thermal expectation value of an operator O is given by 

{d)p = Z-^ tr \d exp(-/3i?)j , (2.2) 

^Scc for example W 



where 

Z = trexp(-/3H) (2.3) 

is the partition function, Tr is the many-body trace, and (3 is the inverse temperature. Ground-state 
(zero-temperature) properties are obtained in the hniit /3 — > cx). 
For a many-body operator H in the quadratic form 

i?-5]6„0„ + i^y„a2 (2.4) 

containing one-body operators Oq. and two-body operators O^ , a HS transformation of the imaginary- 



time evolution operator (2.1) leads to 

1/2 



exp(-/3iJ) « /" n ^^"^" (^^) ^(^) n cxp(- A/3 Ma„)) 



(2.5) 



by splitting /? into Nt time slices such that A/3 = /3/Nt. 

Here, (?„ denotes a set of auxiliary fields (one for each two-body operator appearing in (2.4)) 
at time-slice n, and a stands for the totality of fields at all time-slices. Furthermore, G(cr) is the 
gaussian weight factor 

G(a)=exp(-5]iA/3|K.|0, (2-6) 



a,n 



and h{an) is the one-body Hamiltonian 

h(^n)^'^{^a+ SaVa(7an)Oa (2.7) 

a 

where Sq = ±1 (= ±i) if Vq < (> 0). Thus, the HS transformation has the effect of replacing the 



quadratic dependence on Oa in (2^) with a linear one, at the expense of an integral over auxiliary 
fields. 

Let Ua denote the one-body evolution operator 

f7, =exp(-A/3/i(a„)) (2.8) 

and F[a) its trace: 

F(a)=tr (?/,). (2.9) 

Expectation values can then be written using the above path-integral decomposition of the partition 
function: 

^ _ /P[a]G(a)tr(Q?7.) 



An effective action Sa can be defined such that ( 2.10| ) appears as a simple expectation value: 






with 

^.=5^iA/3|y„K„-logF(a), (2.12) 

where {d)„ = Tv {6 IJ „) / F {a) . 

Further, since exp(— 5o-) is not positive definite, we define the sign-function 



Then, writing e ^'' — ^{<t) e ^■^ with a real Sc, the expectation value (2.10) appears as a ratio of 
expectation values 

^ _ JV[a] e-^- $(a) (O). _ « $.(6). » 

'^ /I?Me-^^$(ff) <*. > ^ ■ ^ 

each calculated with a positive definite probability distribution e^^" . The Monte Carlo average over 
samples is denoted as ^ ... 3>. If the function F{a) is negative on a substantial part of the cr-field 
manifold, the expectation values <C ^a{0)a ^ and <C $,7 ^ can each become very small, and the 
ratio only converges asymptotically. This is the essence of the sign problem in quantum Monte Carlo 
calculations. Below, we construct expectation values susceptible to the sign problem, in order to 
gain insight into the circumstances in which a CL equation approach can be applied successfully. 

3 Complex Langevin Equation and Simple Integrals 

The HS transformation is nothing but the Gaussian integral identity 

e-^"' = ^/ dae-5-%±i--. (3.1) 

V 27r j-oo 



Note that this case corresponds to a repulsive interaction [T^ > in (2.7)]. As (3.1) has no imaginary 
part, this reduces to 



e ^ ^ 



—= daeri''\os{zcT) . (3.2) 

V27r j_oo 



This is precisely the form appearing in the denominator of ( 2.10 ), but in one dimension, with cos(zct) 
playing the role of F{(t). 

Generalizing this, we would like to examine expectation values 

as a function of the real number z and the "particle number" N, using MC evaluation and the 
complex Langevin equation approach. 
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Figure 1: Monte Carlo average of integral (cr^)Ar=i ( |3.3| ) (squares), the exact solution (solid line), 
and the sign of the calculation (insert) 



Clearly, the MC procedure will suffer from the sign-problem only for odd N. Fig. 1 shows a 
straightforward MC evaluation of this integral for A^ = 1, where the inset shows the development 
of the average sign <i>(cr) — cos(z(t)/| cos(z(7)|. As expected, the accuracy of the MC estimate 
deteriorates as ^{a) -^ 0. 

As an alternative to the MC procedure, consider the Langevin equation. For systems with real 
actions, expectation values such as 



with the partition function 



{0)^^JdaO{a)e- 



Z= da e-^(^) 



S(<t) 



(3.4) 



(3.5) 



can be calculated by creating a stochastic process using the Langevin equation, with an equilibrium 
distribution Po(o') = ^~^ exp(— 5(cr)). 
The Langevin equation is given by 



da{t) 



1 dS 



vit), 



(3.6) 



dt 2 da 

where t is a fictitious time and 77 (i) is stochastic noise with zero mean and unit variance: 

{vitMn)=Sit-t'). (3.7) 

To each Langevin equation corresponds a Fokker-Planck equation for the probability density P{a, t) 

dP{a,t) 



dt 



Hfp P{a,t) 



(3.8) 



with a Hermitian Fokker-Planck Hamiltonian 

I d f d dS\ , , 

For solutions with exponential time dependence, 

P(ff,i)=e-^*PB(CT) , (3.10) 

Equation (pjy) reverts to an eigenvalue equation. For i — > oo, P{<7,t) -^ Poi<^), the solution for the 
lowest eigenvalue E — 0. Accordingly, for i ^ oo, if CTo(i) is the solution to the Langevin equation 



(3.e) 



0(f7o(i)) ^ (O) (t^oo) (3.11) 



with (O) given by (3.4). Finally, ergodicity assures that (O) is also obtained by averaging over the 
path ao{t): 



1 '•^ 



(O) = « 0{ao) » - lim - / 0{aoit)) . (3.12) 

In principle, nothing prevents us from using the Langevin equation in the case where the action 
is complex: 

S{a)=SRia)+iSi{c7) . (3.13) 

Then, we obtain two equations, for the real and the imaginary part of a: 

dcTBit) 1 „ fdS\ , , ,„_, 



dt 2 \d(7 



dt 2 \da 

However, the Fokker-Planck Hamiltonian loses its hermeticity, and the eigenvalues can acquire imag- 
inary parts. As a consequence, the probability distribution P{uR,<Ji,t) need not converge anymore, 
nor does the expectation value (O). While the lowest eigenvalue is still Eq = 0, the En with n > 
are in general complex, and the asymptotic condition P{a) —^ Pq is violated whenever there are any 
En with Re En > (n > 0) jl^l. As a rule of thumb then, expectation values obtained via the CL 
equation should only be trusted if the ensemble averages become time-independent pi[ . 



For the numerical solution of equations (3.14,3.15) we use the two-step algorithm of Greenside 



and Helfand |_2[ . Defining the complex gradient as 

^S{t)^ — [aR{t),aj{t)] , (3.16) 

U(T 

the stochastic differential equation is discretized via 

<^Biti/2) = <jR{to)-AtRe[WS{to)] + VAiTj{to) (3.17) 

(TB.{h) = CT7j(io)- 2^*Re[V5(to)+V5(ti/2)] +\/Ai77(io) (3.18) 
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Figure 2: Langevin average of (cr^)jv=i- Squares: extended action, triangles: original action, solid 
line: exact solution. 



and analogously for the imagi nary partQ 

For the expectation value ( |3.3|) , the effective action is 



SN = ^cr^-N log[cos(o-z)] 



while the associated Langevin equation reads 



1 



a = — {a + N z tan(z(T)) + rj 



(3.19) 



(3.20) 



It is immediately apparent that the fixed points of this equation {& = 0) lie on the real axis. On the 



other hand it also clear that, since for example {a 



'N=l 



1 — z , the correct solution (To{t) needs to 



spend a considerable amount of time away from the real line, at least for odd N. Analytically, this 
must happen because (for odd N) there is a delta-function-like drift term that has been ignored in 
(3.20|) which is due to the imaginary part of Sn- The additional drift term is 



SS/Sn — ±i7r(5(cos(zcr)) 



(3.21) 



but cannot be adequately modeled numerically. As a consequence, the solution spends most of its 
time on the real line between the first turning points, and the resulting average is inaccura te, a s has 
been noted previously [Q. The solution to this dilemma is also not new. Going back to (3T) and 
iV = 1 we see that the action can also be written as 



Si 



1 



a ± izcr 



(3.22) 



"^For this particular choice of variables, the corresponding equation for the imaginary part does not include a noise 
term. For other choices, see, e.g., [h3||. 
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Figure 3: Complex Langevin average of (cP')n for N = 2 and N = i. 



in which case the fixed point is away from the real line in the complex plane: ct/ — ±z. Since then 
P{<Tu, c/) — > exp(— i cr|j.) S{(Tj =p z), the complex average reduces to 



{a')N=i 



S^^(j'^P{(JR,(Ti)daRd(Ti JZo^'zz '^^ ^^ exp(-i 0-2) 
/^ ^(o'fl, o-/) dan da I 



n<X}±iz , / 



(3.23) 



The integral over the complex path removed from the real line by ±iz equals the one over the real 
line if the observable has no poles in the enclosed area. 



In Fig. 2 we show the result of a CL evaluation of {(j'^)n=i using the "extended" action (3.22) 



(squares) and the original action ( [3.19| ) (triangles). The solid line represents the exact result. As 
expected, a complex Langevin simulation with fixed points on the real line does not converge, while 
the calculation with the extended action is robust even in the regime where the MC (Fig. 1) fails. 

Can this procedure be extended to arbitrary iV? A canonical extension of the method exists, 
where cos^ {a z) is decomposed into single powers of cosines at multiples of az. Subsequently, the 
cosines are replaced by exponentials. However, this procedure results in a shift of the fixed point 
away from the real line only for odd N . Fig. 3 shows results for the cases N = 2 and iV = 3, with 
exact results 



(ct )n=2 



4^2 



{a')N=3 = l-z' 



l + e2-' ' 



1 



Ip-iz^ 
3^ 



(3.24) 
(3.25) 



For even N the complex Langevin equation does not converge at all z due to the problems described. 
We do not need to worry about this, however, as the MC procedure is very accurate there. 



4 Quantum Mechanical Toy Model 

In this section we show how shifting the fixed points in a CL evaluation of integrals can be used in 
a quantum mechanical model which describes a single shell of angular momentum j with variable 
filling. This toy model is defined by the Hamiltonian 

i? = -ieiv2 + lyj2 (4.1) 

where N is the number operator and J^ the third component of the angular momentum. A HS- 
transformation on the imaginary-time evolution operator yields 

Ttn [e-^") oc /'daodaie-H-o+-?)trjv (eV^*-o-iV^./. ^i^ (42) 

where Ttn is the many-body trace for fixed particle number TV and angular momentum j 

j 
Trjv(O) ~ yj (mi...TOAr|0|TOi...mw) . (4.3) 

'"1 + ---+™JV = -J 

mi ^ ruj 

Since N and Jz commute, the evolution operator does not need to be decomposed into time-slices. 
Also, the contribution from the number operator is a constant factor 

trjv (e-'3ff) = exp(VA^) trjv([/.J (4.4) 

that drops out of the ratios, and will thus be ignored in the following. Above, we defined the 
one-body evolution operator 



L/,^ =e-V/3^J--i . (4.5) 

Defining (f) — \J fi V a\ and as before 

i^Ar(0) = TrAr(C/,J , (4.6) 

we can write the A'^-particle traces in terms of the one-particle trace 

sm((p/2) 
F2(0) = F''{<i>)-F{2<i)) (4.8) 

^3(0) = F\<t>)-?,F{2<t))F{<t>) + 2F{2,4>) (4.9) 

and so on. Note that unlike in the previous section, the even-iV trace is not positive-definite, while 
still being mostly positive. 

We shall focus on the expectation value^ 



/ dcre-^'^'trAr(j2[/^) 



'Note that we now write a instead of cri for simplicity. 



{JDn = . ■ ,i,.p\ , ' (4.10) 



which we rewrite in terms of an effective action as follows 

.1^2 , ^,, -Sjv(<t) 



ai- ^'rry'P-:... . ("d 



where 

and 

1 
2 

For A^ = 1 and arbitrary j, we find for the observable (0 — y/f3Va) 

{Jl)a = J {j + 1) + cot I (u + h cotO' + i)0 - i cot ^ ) (4.14) 



AT^-d^ - log(Fjv(a)). (4.13) 



2 V 2' w ■ 2'^ 2 2 



while the Langevin equation is 



^ = -^^ + ^y^((j + ^)cot(j + i)0-icot0 +7?. (4.15) 

From the oscillatory nature of F((j)) we expect that the MC procedure will become imprecise at 
large /3. In Fig. 4 (left panel) we show the result of a MC calculation of {J'^)p for a. j — 5/2-shell 
and particle numbers N = 1,7V = 2, and TV = 3. Note that TV = 3 corresponds to half-filling, so 
that higher iV's can be described in terms of "hole" -numbers, and revert to the cases displayed. For 
this simple case, the sign (not shown) does not deteriorate too much before the ground state has 
been reached (/3 large), and the calculation is consequently reliable. Let us test nevertheless how a 
complex Langevin approach fares. 

The CL approach based on the Langevin equation ( 4.15| ) with the observable ( 4.14 ) suffers from 



the same problems that we noted with the simple integral: the fixed points are real and the results 
are consequently unreliable. Again, the remedy is to shift the fixed points such that the effective 
path of integration lies in the complex plane. However, here we encounter another difficulty (which 



is also common in more refined shell- model calculations): the expectation value (4.14) has poles (at 
the zeros of F{(j))). Consequently, the expectation value calculated for a shifted path will equal the 
real-path result plus the sum over the poles on the real-line, of which there are infinitely many. For 
this simple toy Hamiltonian, this can be shown to hold true exactly by calculating the residues. In 
more realistic models, however, the sum over the poles is not readily available. Fortunately, for most 
Hamiltonians the observable Trjv (OUa) can be obtained by calculating moments of TrAr Ua- Here, 
for example 

1 d'^ - 

Trw(j2[/^)^__-_^TrArf/, (4.16) 



iSVda^ 



and integrating by parts twice yields 



10 



J = 5/2 
Monte Carlo Longevin 




Figure 4: Average ( J^ ) for a, j — 5/2 shell with N — 1,2,3 obtained with a Monte Carlo approach 
(left panel) and the complex Langevin equation (right panel). For both cases we took 10 samples of 
10,000 updates each. 



With this observable, the action can now be extended into the complex plane. This is achieved as 
in the simple integrals treated previously, by writing the trace in terms of cosines, and replacing 
cos(0) — > exp(i(f>). Quantum mechanically, this amounts to retaining only those terms in the trace 
(4.3) for which the sum of the magnetic quantum numbers mi +m2 + ...m^ (with rrii ^ rrij) is either 



non-negative or non-positive. In this manner, we break time-reversal invariance by hand, since we 
know that it will be taken care of by the symmetry of the integral. Results of the CL calculation are 
shown in the right panel of Fig. 4. The CL averages are stable and accurate at large (3 while accuracy 
deteriorates only for small (3 where large cancellations must occur. Of course, a partial integration 
can also be performed for the Monte Carlo integral, which improves performance markedly because 
of better sampling. In that case, the Monte Carlo results are comparable to those obtained with the 
Langevin equation. 

This changes when the action is forced to acquire substantial complex pieces by "cranking" the 
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Figure 5: Straightforward Monte Carlo average of ( J^ ) for two particles in a cranked j = 5/2 shell 
at /3 = 1, and the sign (inset). 



Hamiltonian (4.1), which corresponds to a toy nucleus with a single shell undergoing a collective 
rotation. Ignoring the term involving the number operator and again setting V2 = 1, the Hamiltonian 
becomes 



He = H - Lo J^ 



(4.18) 



which breaks time-reversal invariance explicitly. Cranking is notoriously difhcult for quantum MC 
calculations because the sign-problem is exacerbated in these cases. In fact, a straightforward MC 
is hopeless because the sign drops very quickly with increasing uj. In Fig. 5 we show ( J^) as well as 
the sign for two particles in a j = 5/2-shell for oj < 1.5, beyond which the MC calculation becomes 
useless. Convergence improves markedly if the observable is simplified by partial integration (Fig. 
6), but of course the sign is still the same and the MC approach fails. 

Le t us calculate this observable with the CL approach. In Fig. 7 we plot the expectation value 
( 4.10| ) as a function of cranking frequency for a, j — 5/2-shell with N = 1,2,3 and (3 — 1.0. While 
the sign in the MC calculation essentially disappears for cj > 3, the accuracy of the CL calculation 
is maintained even as the average sign is small. 



5 Lipkin Model 



To test the CL equation approach in a more realistic situation for a system that can be exactly 
diagonalized, we use the Lipkin model |15 . The Lipkin model is a non-trivial schematic shell model 
capable of describing collective effects in nuclei. It describes N distinguishable particles labeled 
1,2, ...jN, each of which can occupy one of two orbitals with energies (±1/2) (up or down). The 
total number of states therefore is 2^ . 

The Hamiltonian has a one-body term and two two-body terms, and suffers from the sign problem 
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as we shall see below. If written in terms of quasi-spin operators J, it becomes 



H^J,-^v(^Jl + Jl 



Jz-V 



yJx Jyj 



(5.1) 



The eigenstates can be labeled by the total quasi-spin j, and classified into non-degenerate multiplets 
of 2j -I- 1 states each (see, e,g, llB]), from spin j = N/2 down to or i, depending on whether TV is 
even or odd. 

The two-body interaction term does not commute with the free Hamiltonian, which necessitates 
the introduction of time slices in the HS transformation. Writing the imaginary-time evolution 
operator as 



U 



exp(-A/3iJ) 



Nt 



where (3 = A^t A/3, we can apply the Hubbard-Stratonovich transformation to obtain: 
e-^^" ex jva^Vay exp(-iA/3y ^(4")^ + ^(n)2)) n„ exp(-A/3/i^")) , 

■' n 

where 



(5.2) 



(5.3) 



(5.4) 



With this decomposition, exp(— A/3 her) is accurate to order (A/3)^ and U to order A/3. Exact results 
are obtained in the limit A/3 — > 0. 

To obtain averages at finite temperature, we need to take traces of such operators over the many- 
body basis. For small N , the model can easily be diagonalized, which we use to our advantage to 
compare Monte Carlo and Langevin calculations of expectation values with the known exact results. 
We start with the expressions for the static-path approximation (SPA), i..e, for a single time-slice. 
With the representation 



1 

2 

^(o'x +iCTy) 



^ {(^x - icTy) 



V2 



we have 



where 



tri(e-^'"'") = 2 cosh(^ W) , 



W=Jl + 2V^{al~al), 



and indeed quite generally 



Tr7vexp(— A/3ft,o 



2cosh(-^iy) 



N 



(5.5) 



(5.6) 



(5.7) 



(5.8) 
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Similarly, we obtain for example 



(i.)^ = -^tanh(^). (5.9) 



The many-time-slice expressions are straightforward extensions. For time slice n, define 

exp(-A/3/i(,")) = A(") = 4"^ + a(") • r (5.10) 

and 

Ua = exp(-A/34^)), • • • , exp(-A/34^')) = mq + « • r , (5.11) 

where r are the usual Pauli matrices and 

4") = cosh(^T4^„), (5.12) 

a^"^ = ^^sinh(^T4^„) -1^24") | • (5.13) 



V 



Then 



and, for example. 



Trw(;7,) = (2^0)"^ , (5.14) 



J,{a) = TTNiXUa)/TTN{U^) - N^ . (5.15) 

2mo 



The action can be written as 



\ n J 



^- = ^ ( > ] '^i"^' + 4"" 1 - log{TMU,j) , (5.16) 



and observables are obtained as usual. Here, we choose to examine 

The complex Langevin equation requires the gradients of Ti{Ua) with respect to Ux and (jy for 
each time slice: 

5Tr(C/<,) _ ^^ / (1) ^^(„) ^(^^)X ^ ^g^^g^ 



.-, (rt) 

oax 



where Ax is the derivative of the nth slice matrix vl^"). Thus, 

atr(c/. 



4"Ht) = -M4"^-iV^^I#l/tr(C/.) +^ (5.19) 
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Figure 8: Monte Carlo (right panel) and complex Langevin (left panel) average (Jz) in a three 
particle Lipkin model as a function of inverse temperature, with A/3 — 0.05. 



since 



9ai") 



tr(t/^)^ = TVtr (a(i) • • • A'(") • --A^^A ■ Tr(C/<,) 



AT-l 



(5.20) 



if we choose a basis of multi-particle product states, and similarly for Oy. In this manner, all 
remaining traces are over the two-dimensional single-particle space only. Naturally, this equation 
has to be separated into its real and imaginary parts as earlier, since ^X]\i{Uc^ can become negative 
if iV is odd. This is most easily seen examining the single particle trace (p.6[), which becomes 



oscillatory if the argument W = Jl + 2V'^{a1 — Uy) becomes complex. 

In Fig. 8, we show Monte Carlo (right) and complex Langevin (left) calculations of the observable 
{Jz) p as a function of the inverse temperature /3, for the first non-trivial case iV = 3. Even though 
the sign is not strictly positive, the Monte Carlo simulation is very accurate in this case. The 
Langevin averages converge well also, despite the fact that none of the tricks used in the previous 
examples (such as partial integration and extending the action into the complex plane) can be used 
for this model. In order to force a sign problem, we can revert to cranking as before. 

The cranked Hamiltonian 



H — Hn — UJ Jy 



(5.21) 



develops a sign problem because the single particle (and single time-slice) trace (5.6) becomes com- 
plex as the argument (pjh becomes 



W 



2V^{(jl-al)-2^/2Wujay 



1 



(5.22) 



Figure 9 compares Monte Carlo and complex Langevin calculations of the same observable as in 
Fig. 8, at fixed inverse temperature (3 — 2.0 and as a function of the cranking frequency lj. Because 
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Figure 9: (a) Monte Carlo calculation of Re((J2)) as in Fig. 8 (10 samples of 10,000 points), at fixed 
inverse temperate /3 = 2.0 (with A/3 = 0.1) and as a function of cranking frequency w. The inset 
shows the real part of the sign function $i^. (b) Complex Langevin calculation of (J^) with same 
parameters as (a). 



'Yx{Ua) becomes complex (rather than just non-positive), care must be given to the real and imag- 
inary parts of the sign function <&. Indeed, for $ = ^ji -I- i$/ and the observable Oniu) -I- \Oi{p), 
the Monte Carlo average is 



(o> 



/ Vg e-^' {^rOr - $/0/) + i / X'ct e-^' {^jOr + ^rOj) 



r(j). 



(5.23) 



Fig. 9a shows the real part of ( |5.23D as well as the real part of the sign ^r in the inset. The 
sign disappears quickly, both in the MC as well as the Langevin calculation in Fig. 9b, but in the 
Langevin case the accuracy of the average actually increases with increasing cranking frequency. 
Instead, for the MC calculation, this translates into a deteriorating signal-to-noise ratio. Also, the 
complex Langevin calculation does not necessitate the calculation of four separate integrals such as 
in ( 5.23| ). However, for small cranking frequencies oj the complex Langevin averages are noticeably 
non-stationary, which results in larger error bars. 



6 Conclusions 

The complex Langevin equation offers a new perspective on the pervasiveness of the sign problem 
in fermionic quantum many-body calculations. It is not without its own problems, however, most 
notably the absence of a convergence proof of the Langevin averages. The root of non-stationarity 
for some complex Langevin averages lies in the structure of fixed points (attractors) and turning 
points (repellers) in the complex plane. If both the attractors and repellers lie on the real line, it is 
just a matter of time until the trajectory hits a pole in the gradient, and the trajectory is thrown 
far into the complex plane. The first two examples we have treated show how this can be avoided 
by modifying the action such that the fixed points move into the complex plane, without changing 
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the value of the average or the pole structure. The third example, the Lipkin model, showed that 
such a procedure is not necessary if the fixed points are naturally in the complex plane (such as is 
the case at finite cranking frequencies) even though the averages may become non-stationary. In 
these cases, the Langevin equation continues to deliver reliable averages even when the Monte Carlo 
averages have become meaningless. 

While this study certainly suggests that the sign problem can be overcome in particular cases, 
it is by no means certain that the procedure will be as successful in so-called real-life applications, 
with realistic interactions. However, as the payoff is potentially large, we believe that there is now 
enough evidence to try this approach. Another area where this approach deserves to be tested is 
lattice gauge calculations of matter at finite chemical potential /i, which suffer from a sign problem 
because the action becomes complex. As the attractors would naturally move into the complex plane 
at about ^ i/i, the complex Langevin approach seems particularly natural in this case. 

This work was supported in part by the National Science Foundation Grant Nos. PHY9 1-1 5574 
and PHY94-12818, as well as a Cahech Fairchild Fellowship to CA. 
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